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Abstract 

Using chiral perturbation theory, we calculate the density-dependent spin-orbit coupling 
generated by the two-pion exchange three-nucleon interaction involving virtual A-isobar 
excitation. From the corresponding three-loop Hartree and Fock diagrams we obtain an 
isoscalar spin-orbit strength F so (kf) which amounts at nuclear matter saturation density to 
about half of the empirical value of 90 MeVfm 5 . The associated isovector spin-orbit strength 
G so (kf) comes out about a factor of 20 smaller. Interestingly, this three-body spin-orbit 
O , coupling is not a relativistic effect but independent of the nucleon mass M. Furthermore, 

we calculate the three-body spin-orbit coupling generated by two-pion exchange on the 



o 



basis of the most general chiral 7r7T-/ViV-contact interaction. We find similar (numerical) 
results for the isoscalar and isovector spin-orbit strengths F so (kf) and G so (kf) with a 
strong dominance of the p-wave part of the 7T7r./V./V-contact interaction and the Hartree 
contribution. 
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The microscopic understanding the dynamical origin of the strong nuclear spin-orbit force 
is still one of the key problems in nuclear physics. The analogy with the spin-orbit interaction 
in atomic physics gave the hint that it could be a relativistic effect. This idea has lead to the 
construction of the (scalar- vector) mean-field models for nuclear structure calculations [1, 2]. In 
these models the nucleus is described as a collection of independent Dirac quasi-particles moving 
in self-consistently generated scalar and vector mean-fields. The footprints of relativity become 
visible through the large nuclear spin-orbit coupling which emerges in that framework naturally 
from the interplay of the two strong and counteracting (scalar and vector) mean-fields. The 
corresponding many-body calculations are usually carried out in the Hartree approximation, 
ignoring the negative-energy Dirac-sea. The NN-interaction underlying these models is to be 
considered as an effective one that is tailored to properties of finite nuclei but not constrained 
(completely) by the observables of free NN-scattering. 

On the other hand it has long been known that calculations based on Hamiltonians which 
contain only realistic two-nucleon potentials (thus fitting accurately all NN-phase shifts and 
mixing angles below the NN7r-threshold) often cannot predict the observed spin-orbit splittings 
of nuclear levels. In fact one of the original motivations for the Fujita-Miyazawa three-nucleon 
potential [3] was just the study of such spin-orbit splittings. In ref.[4] it has then been shown 
that one out of the Urbana family of three-nucleon forces makes a substantial contribution to 
the spin-orbit splitting in the nucleus 15 N. Moreover, three-nucleon forces are actually needed 
in addition to realistic two-nucleon potentials in order to reproduce the correct saturation point 
of (isospin-symmetric) nuclear matter [5]. The long-range part of the three-nucleon interaction 
is generated in a natural way by two-pion exchange [6] and it can in fact be predicted by using 
chiral symmetry [7]. 

The purpose of this paper is present analytical results for the nuclear spin-orbit coupling 
generated by the (chiral) two-pion exchange three-nucleon interaction. In order to arrive at 
such results we will make use of the density-matrix expansion of Negele and Vautherin [8]. This 
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technique allows one to compute diagrammatically the nuclear energy density functional which 
includes the wanted spin-orbit coupling term proportional to the density gradient. We will first 
consider the three-nucleon interaction proposed originally by Fujita and Miyazawa [3] where 
two pions are exchanged between nucleons while the third nucleon is excited to a (p-wave) A- 
resonance. We are able to evaluate the corresponding three-loop Hartree and Fock diagrams in 
closed analytical form. Then, we will turn to the two-pion exchange three-nucleon interaction 
generated by the (most general) chiral 7T7riViV-contact vertex proportional to the second-order 
low-energy constants Cj. The effects from explicit A-excitation reappear in this description via 
resonance contributions to the low-energy constants 03,4. In both approaches we will separately 
discuss the isoscalar and isovector spin-orbit strengths. 

Let us begin with writing down the explicit form of the spin-orbit coupling term in the nuclear 
energy density functional: 

— * — * — * — * — * — * 

£so[Pp, Pn, J P , Jn] = Vp • J F so (k f ) + Vp v ■ J v G so (k f ) , (1) 

where the sums p = p p + p n , J — J p + J n and differences p v = p p — p n , J v = J p — J n oi proton 
and neutron quantities have been introduced. 

P P ,n(r) = = E ^(0*2(0 , (2) 

" aGocc 

denote the local proton and neutron densities which we have rewritten in terms of the corre- 
sponding (local) proton and neutron Fermi-momenta k PjTl (f) and expressed as sums over the 
occupied single-particle orbitals ^ p a } L {f). The spin-orbit densities of the protons and neutrons 
are defined similarly: 

J P Ar) = £ 9<$(?)i a x W$(f ) . (3) 

agocc 

Furthermore, F so (kf) and G so (kf) in eq.(l) denote the density dependent isoscalar and isovector 
spin-orbit strength functions. In Skyrme parameterizations [9] these are just constants, F so (kf) = 
3G so (kf) = 3W /4, whereas in our calculation their explicit density dependence originates from 
the finite range character of the two-pion exchange three-nucleon interaction. 

The starting point for the construction of an explicit nuclear energy density functional 
£ so [. . .] is the bilocal density- matrix as given by a sum over the occupied energy eigenfunc- 
tions: J2aeocc ^p°n(f ~~ Q-/2)^ p a ^ (r + a/2). According to Negele and Vautherin [8] it can be 
expanded in relative and center-of-mass coordinates, a and r, with expansion coefficients deter- 
mined by purely local quantities (nucleon density, kinetic energy density and spin-orbit density). 
As outlined in Sec. 2 of ref. [10] the Fourier-transform of the (so expanded) density-matrix defines 
in momentum-space a medium-insertion T(p,q) for the inhomogeneous many-nucleon system. 
It is straightforward to generalize this construction to the isospin-asymmetric situation of dif- 
ferent proton and neutron local densities p p<n {f ) an d J p , n (f). We display here only that part of 
the medium-insertion T(p,q) which is actually relevant for the diagrammatic calculation of the 
isoscalar and isovector spin-orbit terms defined in eq.(l): 

T(p,q) = J dhe~^[ ] ^e{k p -\p\)+ l -^e{k n -\p\) 

'S(k f - \p\) - k f 6'(k f - \p\)] (a x p) ■ (J + r 3 J v )] . (4) 
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When working to quadratic order in deviations from isospin symmetry (i.e. proton-neutron 
differences) it is sufficient to use an average Fermi-momentum kf in the prefactor of the spin-orbit 

— * — * 

density J + r 3 J v . The double-dash in the left picture of Fig. 1 symbolizes the medium insertion 



T(p,q) together with the assignment of the out- and in-going nucleon momenta p ± q/2. The 
momentum transfer q is provided by the Fourier-components of the inhomogeneous (matter) 
distributions p Pi „(r) and J Pi „(r ). 




Fig. 1: Left: The double-dash symbolizes the medium insertion T(p,q) defined by eq.(4)- Next 
shown are the three-loop two-pion exchange Hartree and Fock diagrams involving one chiral 
nnNN -contact vertex (symbolized by the heavy dot). The combinatoric factors of these diagrams 
are 1/2 and 1, in the order shown. 




Fig. 2: Two-pion exchange Hartree and Fock diagrams with (single) virtual A-isobar excitation. 
The solid double-line denotes the A-isobar and dashed and solid lines represent pions and nu- 
cleons, respectively. For is o spin- symmetric nuclear matter the isospin factors of these diagrams 
are 8, 0, and 8. The combinatoric factor is 1 in each case. 

Now we turn to the analytical evaluation of the two-pion exchange diagrams with (single) 
A-isobar excitations shown in Fig. 2. We give for each diagram only the final result for the spin- 
orbit strengths F so (kf) and G so (kf) omitting all technical details related to extensive algebraic 
manipulations and solving elementary integrals. Putting a medium insertion at each of the three 
nucleon propagators of the Hartree diagram (left diagram in Fig. 2) we obtain the following 
contribution to the isoscalar spin-orbit strength: 

so[f) 8tt 2 A/4 1 1 + 4u* 4u [ +4M J/' W 

where u = kf/m n denotes the ratio of the two small scales kf and m^. The A-propagator shows 
up in this expression merely via the (reciprocal) AiV-mass splitting, A = 293 MeV. Additional 
corrections to the A-propagator coming from differences of nucleon kinetic energies etc. will make 
a contribution to the spin-orbit strength F so (kf) at least one order higher in the small momentum 
expansion. In eq.(5) we have inserted the empirically well-satisfied relation g^NA = 



for the 7riVA-coupling constant together with the Goldberger-Treiman relation = g^M/ f v . 
Let us briefly sketch the main mechanism which generates the strength function F so (kf). The 
exchanged pion-pair in the Hartree digram transfers a momentum q between the left and the 
right nucleon ring and this momentum q enters also the pseudovector 7riViV-interaction vertices. 
The spin-orbit strength F so (kf) arises from the spin-trace tr[<j- (Q + q/2) a- (Q — q/2) a-(pxj)] = 
2i(q x Q) ■ (p x J) where iq gets converted to V/c/ = (n 2 / 2kj)V p by Fourier transformation. 
The rest is a solvable integral over the product of three Fermi spheres. The second Fock diagram 
in Fig. 2 (with parallel pion lines) has the isospin factor for isospin-symmetric nuclear matter 
and from the third Fock diagram (with crossed pion lines) we get the following contribution to 
the isoscalar spin-orbit strength: 

F so (k f Y A - Fock) = S^2y I 8 " 2 " 12 + (3m " + 4) ln(1 + 4 " 2) ' 2 • (6) 

It is highly remarkable that the pertinent nine-dimensional integral over the product of two 
(different) pion-propagators and other momentum dependent factors can solved in terms of (a 
square of) elementary functions without the occurrence of any dilogarithm. The specific isospin 
structures of the irNN- and 7riVA-vertices determine uniquely the ratio of isovector to isoscalar 
spin-orbit strength of each of the three diagrams. We find that the (left) Hartree diagram in 
Fig. 2 does not contribute to the isovector spin-orbit strength G so (kf) while the combined result 
of both Fock diagrams in Fig. 2 reads: 

G S o(M (A " Fock) = ^so(M (A " F ° ck) , (7) 

with a contribution of the second and third Fock diagram in the ratio 6:1. It is important to 
note that the expressions in eqs.(5,6) are independent of the nucleon mass M and therefore these 
27r-exchange three-body spin-orbit couplings are not relativistic effects. In fact the diagrams in 
Fig. 2 with two medium insertions (on non-neighboring nucleon propagators) do also generate 
two-body spin-orbit couplings. The latter are however genuine relativistic effects proportional 
to 1/M and therefore counted as one order higher in the small momentum expansion. Such two- 
body contributions to the spin-orbit strengths F so (kf) and G so (kf) can generally be expressed 
in terms of the spin-orbit amplitudes entering the T-matrix of elastic NN-scattering: 

F so (k f ) {2 ~ hody) = ~{3Vso(0) + Vso(2k f ) + 3W S o(2k f ) 

+ dxx[v so (2xk f ) +3W so (2xk f )]Y (8) 

Gso(A;/) (2_body) = l{w S o(2k f )-Vso(2k f )-3Wso(0) 

+ jf 1 dxx[w so (2xk f ) - V S o(2xkf)] | . (9) 

The terms — V S o(0)/2 and — W S o(0)/2 belong to Hartree- type diagrams (with two closed nucleon 
lines) while the remaining ones summarize the contributions from Fock-type diagrams (having 
just one closed nucleon line). Explicit expressions for the isoscalar and isovector spin-orbit 
NN-amplitudes Vso{q) and Wso{q) as they arise from 27r-exchange with (single and double) 
A-excitation can be found in the appendix of ref.fll] (modulo regularization dependent additive 
constants). 

For the numerical evaluation of eqs. (5,6,7) we use the (physical) parameters: M = 939 MeV 
(nucleon mass), = 135 MeV (neutral pion mass), f w = 92.4 MeV (pion decay constant) and 
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qa = 1-3 (equivalent to a 7riViV-coupling constant of g n N = 9aM/ f n = 13.2). The full line 
in Fig. 3 shows the isoscalar spin-orbit strength F so (kf) generated by the two-pion exchange 
three-nucleon interaction involving virtual A-excitation as a function of the nucleon density 
p = 2/cj/37r 2 . As it is typical for a three-body effect the spin-orbit strength F so (kf) starts from 
the value zero at zero density p = 0. The contribution of the Hartree diagram is by far the 
dominant one. At nuclear matter saturation density (where kf ~ 2m„-) one finds for example 
F so (2m^)( A - Hart ) = 48.2 MeVfm 5 to be compared with a Fock contribution of F so (2m u )( A - Fock ) = 
1.2MeVfm 5 . Clearly, this 27r-exchange three-body spin-orbit coupling is sizeable [4]. In the 
region around saturation density p ^ 0.17 fm 3 it amounts to about half of the "empirical" 
value 3W /4 ~ 90MeVfm 5 deduced in the Skyrme phenomenology of nuclear structure [9]. 
The findings of refs.[4, 6] concerning spin-orbit splittings in light nuclei point of course in the 
same direction. Finally, the dashed line in Fig. 3 shows the isovector spin-orbit strength G so (kf) 
(magnified by a factor 10). In comparison to the isoscalar spin-orbit strength F so (kf) it is only 
a small 5% correction. 
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Fig. 3: The spin-orbit strength generated by the two-pion exchange three-nucleon interaction 
involving virtual A-isobar excitation versus the nucleon density p = 2kj/3n 2 . The full 



curve 



shows the isoscalar spin-orbit strength F so (kf) and the dashed curve shows the isovector spin- 
orbit strength G so (kf) magnified by a factor 10. 

Next, we turn to a more general derivation of the three-body spin-orbit coupling generated 
by two-pion exchange. Chiral symmetry determines the 27r-exchange three-nucleon interaction 
uniquely at leading order [7] . It follows from a tree-diagram involving the chiral 7T7rA^A^-contact 
vertex proportional to the second-order low-energy constants cf 



■pi 2cU 

J TV 



c 3 q a -q b - 2ci ml + c 4 e abc r c a ■ (q a x q b 



(10) 



Here, q a ^ denote out-going pion-momenta and we have already dropped the C2-term proportional 
to the product of two pion-energies. In the present application these (off-shell) pion-energies 
become equal to differences of nucleon kinetic energies thus producing a higher order relativistic 
l/M 2 -correction. A straightforward evaluation of the (three-loop) Hartree diagram in Fig. 1 with 
three medium insertions gives the following contribution to the isoscalar spin-orbit strength: 



F so (kf) 



(cj— Hart) 



47r2/ 4 



( Cl - C 3) W -2 C3 ^ + C3-C lln(1+4u2) 



1 +4m 2 



4m 



where u = kf/m^. The contributions of the s-wave part and the p-wave part of the irirNN- 
contact vertex are distinguished by the two (isoscalar) low-energy constants c\ and C3. The 
isovectorial and spin-dependent C4-vertex makes a non-vanishing contribution only through the 
Fock diagram. We find the following total result for the isoscalar spin-orbit strength from the 
three-loop Fock diagram in Fig. 1: 



^so(fc/) 



(cj— Fock) 



4vr 2 (8A) 4 
-(c 3 + c 4 ) 



32ci 



4m 



u 



-1 



ln(l + 4m 2 ) 



8u 2 



12 + (3m -2 + 4) ln(l + 4m 2 



(12) 



Its analytical form is remarkably simple. Again there no contribution from the Hartree diagram 
in Fig. 1 to the isovector spin-orbit strength G so (kf) while the Fock diagram in Fig. 1 leads to 
the combined result: 



G so (k f ) 



(cj —Fock) 



47T 2 (8/ 7r ) 4 



1 2 



-32ci 4u-iT 1 m(l + Au 2 ) 



8u 2 - 12 + (3m~ 2 + 4) ln(l + 4m 2 ) 



(13) 



Here, the isospin structures of the ci^-vertex and the Czt-vertex shows up through relative fac- 
tors 1 and —1/3. We note aside that the previously calculated contributions from explicit 
A-excitation (see eqs. (5,6,7)) reappear in the (general) expressions eqs. (11,12,13) through A- 
resonance contributions to the low-energy constants: — c 3 A ^ = 2c^ = g\/2A. This connection 
allows one also to trace back the origin of the factor 7/3 in eq.(7). Finally, we consider the 
(leading-order) Weinberg- Tomozawa 7r7riViV-contact vertex. It generates a three-body spin-orbit 
coupling first in form of a relativistic 1/M-correction. Because of the isovector nature of the 
Weinberg- Tomozawa contact- vertex the Hartree diagram (in Fig. 1) can make a contribution 
only to the isovector spin-orbit strength: 



G so (k f ) 



(WT-Hart) 



3 arctan 2u — 2u — u 1 ln(l + 4m 2 ) 



(14) 



96tt 2 M/4 

The Fock diagram (in Fig. 1) on the other hand generates isoscalar and isovector spin-orbit 
strengths in a fixed (relative) ratio: 

F so (A; / )( WT - Fock ) = -3G so (A; / )( WT - Fock ) 



g 2 A m w u 2 
" A^M{AUf y 
3 + 2u 2 + 10m 4 
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15m 2 arctan 2m H u — 39m 3 

4m 

, /-, . 2\ 3 + 8m 2 -64m 4 , 2 2 , 
In 1 + 4m 2 + -— In 2 1 + 4m 2 

v ' 64m 5 



3(m 2 + m 4 ) + 6(m + m 3 )(6x 2 - 1 - m 2 )L 



3(1 + m ) - (15 + 34m 2 + 19m 4 )x 2 + (33 + 29m 2 )x 4 - 13x fc 



(15) 



with the auxiliary function: 

4x 1 + (u - x) 2 y ' 

The power of the pion mass m n in their prefactors indicates that all contributions written in 
eqs. (11-15) are of the same order in the small momentum expansion. 
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Fig. 4-' The spin-orbit strength generated by the two-pion exchange three-nucleon interaction 
involving the general chiral titiNN -contact vertex versus the nucleon density p = 2kj/3n 2 . The 
full curve shows the isoscalar spin-orbit strength F so (kf) and the dashed curve shows the isovector 
spin-orbit strength G so (kf) magnified by a factor 10. 

For the numerical evaluation we use the values c\ = — 0.64GeV _1 , C3 = — 3.90GeV _1 and 
C4 = 2.25 GeV -1 of the low-energy constants which have been determined (at tree-level) in 
ref.[12] from some low-energy 7riV-data. The full line in Fig. 4 shows the resulting isoscalar 
spin-orbit strength F so (kf) as a function of the nucleon density p = 2kj/3n 2 . The contribution 
of the Hartree diagram involving the p-wave contact vertex proportional to c 3 is the absolutely 
dominant one. For example one finds at nuclear matter saturation density (where kf ~ 2m n ) the 
value F so (2m 7r )( C3 - Hart ) = 65.3MeVfm 5 . In comparison to this the s-wave Hartree contribution is 
very small, F so (2m 7r ) ( ^ Cl ~ Hart ^ = 3.6MeVfm 5 . Moreover, the even smaller Fock contributions have 
a tendency of cancelling each other. The somewhat smaller values of F so (kf) in Fig. 3 compared 
to those in Fig. 4 originate mainly from the fact that the A-resonance saturates the low-energy 
constant C3 to about only three quarters in magnitude: c^ = —g\/2A ~ — 2.9GeV _1 . The 
dashed line in Fig. 4 shows the isovector spin-orbit strength G so (kf) (magnified by a factor 10) as 
a function of the nucleon density p = 2kj/3ir 2 . Again, in comparison to the isoscalar spin-orbit 
strength F so (kf) it is only a small 5% correction. The largest individual contribution comes here 



from the Fock diagram involving the p-wave C3-contact vertex which for example gives at nuclear 
matter saturation density: G , so (2m 7r )( C3_Fock ) = 3.2MeVfm 5 . 

In summary we have calculated in this work the spin-orbit coupling generated by the two- 
pion exchange three-nucleon interaction. We have made use of the density-matrix expansion of 
Negele and Vautherin [8]. This method allows one to compute diagrammatically the nuclear 
energy density functional which contains the spin-orbit coupling term of interest. We have 
derived simple analytical expressions for the density- dependent isoscalar and isovector spin-orbit 
strengths F so (kf) and G so (kf). First, we have considered the two-pion exchange three-nucleon 
interaction of Fujita and Miyazawa [3] where one nucleon is excited to a p-wave A- resonance. 
The corresponding three-loop Hartree and Fock diagrams generate spin-orbit couplings which 
are not relativistic effects but independent of the nucleon mass. The Hartree diagram and 
the isoscalar component F so (kf) are by far dominant. At nuclear matter saturation density 
these A-driven three-body mechanisms generate about half of the empirical isoscalar spin-orbit 
strength. The calculations of spin-orbit splittings in light nuclei [4, 6] point of course in the same 
direction. Secondly, have we derived more generally the three-body spin-orbit coupling generated 
by two-pion exchange on the basis of the chiral 7T7rAA-contact vertex. In that framework we 
have obtained similar (numerical) results for the density dependent isoscalar and isovector spin- 
orbit strengths F so (kf) and G so (kf). The p-wave part of the chiral 7T7rAA-contact interaction 
(proportional to the low-energy constant C3) and the Hartree diagram give rise to the absolutely 
dominant contribution. 

On the other hand it has been shown recently in ref. [10] that iterated one-pion exchange 
generates an isoscalar spin-orbit strength F so (kf) that is sizeable but of the wrong negative sign. 
Combining those results [10] with the present ones, one may conclude that the net nuclear spin- 
orbit coupling generated by (multi) pion-exchange is rather small, at least for densities around 
nuclear matter saturation density. Lorentz scalar and vector mean-fields with their in-medium 
behavior governed by QCD sum rules could therefore be the appropriate dynamical framework 
for building up the strong (isoscalar) nuclear spin-orbit interaction. Indeed such a proposal has 
recently been successfully applied in ref. [13] to nuclear structure calculations. 
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